Fossil-calibrated molecular clock data enable reconstruction of steps leading to differentiated multicellularity and anisogamy in the Volvocine algae.

BACKGROUND: Throughout its nearly four-billion-year history, life has undergone evolutionary transitions in which simpler subunits have become integrated to form a more complex whole. Many of these transitions opened the door to innovations that resulted in increased biodiversity and/or organismal efficiency. The evolution of multicellularity from unicellular forms represents one such transition, one that paved the way for cellular differentiation, including differentiation of male and female gametes. A useful model for studying the evolution of multicellularity and cellular differentiation is the volvocine algae, a clade of freshwater green algae whose members range from unicellular to colonial, from undifferentiated to completely differentiated, and whose gamete types can be isogamous, anisogamous, or oogamous. To better understand how multicellularity, differentiation, and gametes evolved in this group, we used comparative genomics and fossil data to establish a geologically calibrated roadmap of when these innovations occurred. RESULTS: Our ancestral-state reconstructions, show that multicellularity arose independently twice in the volvocine algae. Our chronograms indicate multicellularity evolved during the Carboniferous-Triassic periods in Goniaceae + Volvocaceae, and possibly as early as the Cretaceous in Tetrabaenaceae. Using divergence time estimates we inferred when, and in what order, specific developmental changes occurred that led to differentiated multicellularity and oogamy. We find that in the volvocine algae the temporal sequence of developmental changes leading to differentiated multicellularity is much as proposed by David Kirk, and that multicellularity is correlated with the acquisition of anisogamy and oogamy. Lastly, morphological, molecular, and divergence time data suggest the possibility of cryptic species in Tetrabaenaceae. CONCLUSIONS: Large molecular datasets and robust phylogenetic methods are bringing the evolutionary history of the volvocine algae more sharply into focus. Mounting evidence suggests that extant species in this group are the result of two independent origins of multicellularity and multiple independent origins of cell differentiation. Also, the origin of the Tetrabaenaceae-Goniaceae-Volvocaceae clade may be much older than previously thought. Finally, the possibility of cryptic species in the Tetrabaenaceae provides an exciting opportunity to study the recent divergence of lineages adapted to live in very different thermal environments.

Different transcription of novel, functional long non-coding RNA genes by UV-B in green algae, Volvox carteri.

Long non-coding RNAs (lncRNAs) are identified as important regulatory molecules related to diverse biological processes. In recent years, benefiting from the rapid development of high-throughput sequencing technology, RNA-seq, and analysis methods, more lncRNAs have been identified and discovered in various plant and algal species. However, so far, only limited studies related to algal lncRNAs are available. Volvox carteri f. nagariensis is the best multicellular model organism to study in developmental and evolutionary biology; therefore, studying and increasing information about this species is important. This study identified lncRNAs in the multicellular green algae Volvox carteri and 1457 lncRNAs were reported, using RNA-seq data and with the help of bioinformatics tools and software. This study investigated the effect of low-dose UV-B radiation on changes in the expression profile of lncRNAs in gonidial and somatic cells. The differential expression of lncRNAs was analyzed between the treatment (UV-B) and the control (WL) groups in gonidial and somatic cells. A total of 37 and 26 lncRNAs with significant differential expression in gonidial and somatic cells, respectively, were reported. Co-expression analysis between the lncRNAs and their neighbor protein-coding genes (in the interval of ± 10 Kb) was accomplished. In gonidial cells, 184 genes with a positive correlation and 13 genes with a negative correlation (greater than 0.95), and in somatic cells, 174 genes with a positive correlation, and 18 genes with a negative correlation were detected. Functional analysis of neighboring coding genes was also performed based on gene ontology. The results of the current work may help gain deeper insight into the regulation of gene expression in the studied model organism, Volvox carteri.

Identification of cell-type specific alternative transcripts in the multicellular alga Volvox carteri.

BACKGROUND: Cell type specialization is a hallmark of complex multicellular organisms and is usually established through implementation of cell-type-specific gene expression programs. The multicellular green alga Volvox carteri has just two cell types, germ and soma, that have previously been shown to have very different transcriptome compositions which match their specialized roles. Here we interrogated another potential mechanism for differentiation in V. carteri, cell type specific alternative transcript isoforms (CTSAI). METHODS: We used pre-existing predictions of alternative transcripts and de novo transcript assembly with HISAT2 and Ballgown software to compile a list of loci with two or more transcript isoforms, identified a small subset that were candidates for CTSAI, and manually curated this subset of genes to remove false positives. We experimentally verified three candidates using semi-quantitative RT-PCR to assess relative isoform abundance in each cell type. RESULTS: Of the 1978 loci with two or more predicted transcript isoforms 67 of these also showed cell type isoform expression biases. After curation 15 strong candidates for CTSAI were identified, three of which were experimentally verified, and their predicted gene product functions were evaluated in light of potential cell type specific roles. A comparison of genes with predicted alternative splicing from Chlamydomonas reinhardtii, a unicellular relative of V. carteri, identified little overlap between ortholog pairs with alternative splicing in both species. Finally, we interrogated cell type expression patterns of 126 V. carteri predicted RNA binding protein (RBP) encoding genes and found 40 that showed either somatic or germ cell expression bias. These RBPs are potential mediators of CTSAI in V. carteri and suggest possible pre-adaptation for cell type specific RNA processing and a potential path for generating CTSAI in the early ancestors of metazoans and plants. CONCLUSIONS: We predicted numerous instances of alternative transcript isoforms in Volvox, only a small subset of which showed cell type specific isoform expression bias. However, the validated examples of CTSAI supported existing hypotheses about cell type specialization in V. carteri, and also suggested new hypotheses about mechanisms of functional specialization for their gene products. Our data imply that CTSAI operates as a minor but important component of V. carteri cellular differentiation and could be used as a model for how alternative isoforms emerge and co-evolve with cell type specialization.

Cell Type-Specific Promoters of Volvox carteri for Molecular Cell Biology Studies.

The multicellular green alga Volvox carteri has emerged as a valuable model organism for investigating various aspects of multicellularity and cellular differentiation, photoreception and phototaxis, cell division, biogenesis of the extracellular matrix and morphogenetic movements. While a range of molecular tools and bioinformatics resources have been made available for exploring these topics, the establishment of cell type-specific promoters in V. carteri has not been achieved so far. Therefore, here, we conducted a thorough screening of transcriptome data from RNA sequencing analyses of V. carteri in order to identify potential cell type-specific promoters. Eventually, we chose two putative strong and cell type-specific promoters, with one exhibiting specific expression in reproductive cells (gonidia), the PCY1 promoter, and the other in somatic cells, the PFP promoter. After cloning both promoter regions, they were introduced upstream of a luciferase reporter gene. By using particle bombardment, the DNA constructs were stably integrated into the genome of V. carteri. The results of the expression analyses, which were conducted at both the transcript and protein levels, demonstrated that the two promoters drive cell type-specific expression in their respective target cell types. Transformants with considerably diverse expression levels of the chimeric genes were identifiable. In conclusion, the screening and analysis of transcriptome data from RNA sequencing allowed for the identification of potential cell type-specific promoters in V. carteri. Reporter gene constructs demonstrated the actual usability of two promoters. The investigated PCY1 and PFP promoters were proven to be potent molecular tools for genetic engineering in V. carteri.

The Genetics of Fitness Reorganization during the Transition to Multicellularity: The Volvocine regA-like Family as a Model.

The evolutionary transition from single-celled to multicellular individuality requires organismal fitness to shift from the cell level to a cell group. This reorganization of fitness occurs by re-allocating the two components of fitness, survival and reproduction, between two specialized cell types in the multicellular group: soma and germ, respectively. How does the genetic basis for such fitness reorganization evolve? One possible mechanism is the co-option of life history genes present in the unicellular ancestors of a multicellular lineage. For instance, single-celled organisms must regulate their investment in survival and reproduction in response to environmental changes, particularly decreasing reproduction to ensure survival under stress. Such stress response life history genes can provide the genetic basis for the evolution of cellular differentiation in multicellular lineages. The regA-like gene family in the volvocine green algal lineage provides an excellent model system to study how this co-option can occur. We discuss the origin and evolution of the volvocine regA-like gene family, including regA-the gene that controls somatic cell development in the model organism Volvox carteri. We hypothesize that the co-option of life history trade-off genes is a general mechanism involved in the transition to multicellular individuality, making volvocine algae and the regA-like family a useful template for similar investigations in other lineages.

A personal cost of cheating can stabilize reproductive altruism during the early evolution of clonal multicellularity.

Understanding how cooperation evolved and is maintained remains an important and often controversial topic because cheaters that reap the benefits of cooperation without paying the costs can threaten the evolutionary stability of cooperative traits. Cooperation-and especially reproductive altruism-is particularly relevant to the evolution of multicellularity, as somatic cells give up their reproductive potential in order to contribute to the fitness of the newly emerged multicellular individual. Here, we investigated cheating in a simple multicellular species-the green alga Volvox carteri, in the context of the mechanisms that can stabilize reproductive altruism during the early evolution of clonal multicellularity. We found that the benefits cheater mutants can gain in terms of their own reproduction are pre-empted by a cost in survival due to increased sensitivity to stress. This personal cost of cheating reflects the antagonistic pleiotropic effects that the gene coding for reproductive altruism-regA-has at the cell level. Specifically, the expression of regA in somatic cells results in the suppression of their reproduction potential but also confers them with increased resistance to stress. Since regA evolved from a life-history trade-off gene, we suggest that co-opting trade-off genes into cooperative traits can provide a built-in safety system against cheaters in other clonal multicellular lineages.

Cellular organization in lab-evolved and extant multicellular species obeys a maximum entropy law.

The prevalence of multicellular organisms is due in part to their ability to form complex structures. How cells pack in these structures is a fundamental biophysical issue, underlying their functional properties. However, much remains unknown about how cell packing geometries arise, and how they are affected by random noise during growth - especially absent developmental programs. Here, we quantify the statistics of cellular neighborhoods of two different multicellular eukaryotes: lab-evolved 'snowflake' yeast and the green alga Volvox carteri. We find that despite large differences in cellular organization, the free space associated with individual cells in both organisms closely fits a modified gamma distribution, consistent with maximum entropy predictions originally developed for granular materials. This 'entropic' cellular packing ensures a degree of predictability despite noise, facilitating parent-offspring fidelity even in the absence of developmental regulation. Together with simulations of diverse growth morphologies, these results suggest that gamma-distributed cell neighborhood sizes are a general feature of multicellularity, arising from conserved statistics of cellular packing.

Molecular and cellular dynamics of early embryonic cell divisions in Volvox carteri.

Cell division is fundamental to all organisms and the green alga used here exhibits both key animal and plant functions. Specifically, we analyzed the molecular and cellular dynamics of early embryonic divisions of the multicellular green alga Volvox carteri (Chlamydomonadales). Relevant proteins related to mitosis and cytokinesis were identified in silico, the corresponding genes were cloned, fused to yfp, and stably expressed in Volvox, and the tagged proteins were studied by live-cell imaging. We reveal rearrangements of the microtubule cytoskeleton during centrosome separation, spindle formation, establishment of the phycoplast, and generation of previously unknown structures. The centrosomes participate in initiation of spindle formation and determination of spindle orientation. Although the nuclear envelope does not break down during early mitosis, intermixing of cytoplasm and nucleoplasm results in loss of nuclear identity. Finally, we present a model for mitosis in Volvox. Our study reveals enormous dynamics, clarifies spatio-temporal relationships of subcellular structures, and provides insight into the evolution of cell division.

Cell Type-Specific Pherophorins of Volvox carteri Reveal Interplay of Both Cell Types in ECM Biosynthesis.

The spheroidal green algae Volvox carteri serves as a model system to investigate the formation of a complex, multifunctional extracellular matrix (ECM) in a relatively simple, multicellular organism with cell differentiation. The V. carteri ECM is mainly composed of hydroxyproline-rich glycoproteins (HRGPs) and there are diverse region-specific, anatomically distinct structures in the ECM. One large protein family with importance for ECM biosynthesis stands out: the pherophorins. The few pherophorins previously extracted from the ECM and characterized, were specifically expressed by somatic cells. However, the localization and function of most pherophorins is unknown. Here, we provide a phylogenetic analysis of 153 pherophorins of V. carteri and its unicellular relative Chlamydomonas reinhardtii. Our analysis of cell type-specific mRNA expression of pherophorins in V. carteri revealed that, contrary to previous assumptions, only about half (52%) of the 102 investigated pherophorin-related genes show stronger expression in somatic cells, whereas about one-third (34%) of the genes show significant higher expression in reproductive cells (gonidia). We fused two pherophorin genes that are expressed by different cell types to yfp, stably expressed them in Volvox and studied the tagged proteins by live-cell imaging. In contrast to earlier biochemical approaches, this genetic approach also allows the in vivo analysis of non-extractable, covalently cross-linked ECM proteins. We demonstrate that the soma-specific pherophorin SSG185 is localized in the outermost ECM structures of the spheroid, the boundary zone and at the flagellar hillocks. SSG185:YFP is detectable as early as 1.5 h after completion of embryogenesis. It is then present for the rest of the life cycle. The gonidia-specific pherophorin PhG is localized in the gonidial cellular zone 1 ("gonidial vesicle") suggesting its involvement in the protection of gonidia and developing embryos until hatching. Even if somatic cells produce the main portion of the ECM of the spheroids, ECM components produced by gonidia are also required to cooperatively assemble the total ECM. Our results provide insights into the evolution of the pherophorin protein family and convey a more detailed picture of Volvox ECM synthesis.

Genome sequencing of the multicellular alga Astrephomene provides insights into convergent evolution of germ-soma differentiation.

Germ-soma differentiation evolved independently in many eukaryotic lineages and contributed to complex multicellular organizations. However, the molecular genetic bases of such convergent evolution remain unresolved. Two multicellular volvocine green algae, Volvox and Astrephomene, exhibit convergent evolution of germ-soma differentiation. The complete genome sequence is now available for Volvox, while genome information is scarce for Astrephomene. Here, we generated the de novo whole genome sequence of Astrephomene gubernaculifera and conducted RNA-seq analysis of isolated somatic and reproductive cells. In Volvox, tandem duplication and neofunctionalization of the ancestral transcription factor gene (RLS1/rlsD) might have led to the evolution of regA, the master regulator for Volvox germ-soma differentiation. However, our genome data demonstrated that Astrephomene has not undergone tandem duplication of the RLS1/rlsD homolog or acquisition of a regA-like gene. Our RNA-seq analysis revealed the downregulation of photosynthetic and anabolic gene expression in Astrephomene somatic cells, as in Volvox. Among genes with high expression in somatic cells of Astrephomene, we identified three genes encoding putative transcription factors, which may regulate somatic cell differentiation. Thus, the convergent evolution of germ-soma differentiation in the volvocine algae may have occurred by the acquisition of different regulatory circuits that generate a similar division of labor.

Characterization and Transformation of reg Cluster Genes in Volvox powersii Enable Investigation of Convergent Evolution of Cellular Differentiation in Volvox.

The evolution of germ-soma cellular differentiation represents a key step in the evolution of multicellular individuality. Volvox carteri and its relatives, the volvocine green algae, provide a model system for studying the evolution of cellular differentiation. In V. carteri, the regA gene controls somatic cell differentiation and is found in a group of paralogs called the reg cluster, along with rlsA, rlsB, and rlsC. However, the developmental program of V. carteri is derived compared to other volvocine algae. Here we examine Volvox powersii which possesses an ancestral developmental program and independent evolution of the Volvox body plan. We sequenced the reg cluster from V. powersii wild-type and a mutant with fewer cells and altered germ-soma ratio. We found that the mutant strain's rlsB gene has a deletion predicted to cause a truncated protein product. We developed a genetic transformation procedure to insert wild-type rlsB into the mutant strain. Transformation did not result in phenotypic rescue, suggesting the rlsB mutation is insufficient for generating the mutant phenotype. The transformation techniques and sequences described here provide essential tools to study V. powersii, a species well suited for studying the evolution of cellular differentiation and convergent evolution of Volvox morphology.

Phylotranscriptomics points to multiple independent origins of multicellularity and cellular differentiation in the volvocine algae.

BACKGROUND: The volvocine algae, which include the single-celled species Chlamydomonas reinhardtii and the colonial species Volvox carteri, serve as a model in which to study the evolution of multicellularity and cellular differentiation. Studies reconstructing the history of this group have by and large relied on datasets of one to a few genes for phylogenetic inference and ancestral character state reconstruction. As a result, volvocine phylogenies lack concordance depending on the number and/or type of genes (i.e., chloroplast vs nuclear) chosen for phylogenetic inference. While multiple studies suggest that multicellularity evolved only once in the volvocine algae, that each of its three colonial families is monophyletic, and that there have been at least three independent origins of cellular differentiation in the group, other studies call into question one or more of these conclusions. An accurate assessment of the evolutionary history of the volvocine algae requires inference of a more robust phylogeny. RESULTS: We performed RNA sequencing (RNA-seq) on 55 strains representing 47 volvocine algal species and obtained similar data from curated databases on 13 additional strains. We then compiled a dataset consisting of transcripts for 40 single-copy, protein-coding, nuclear genes and subjected the predicted amino acid sequences of these genes to maximum likelihood, Bayesian inference, and coalescent-based analyses. These analyses show that multicellularity independently evolved at least twice in the volvocine algae and that the colonial family Goniaceae is not monophyletic. Our data further indicate that cellular differentiation arose independently at least four, and possibly as many as six times, within the volvocine algae. CONCLUSIONS: Altogether, our results demonstrate that multicellularity and cellular differentiation are evolutionarily labile in the volvocine algae, affirming the importance of this group as a model system for the study of major transitions in the history of life.

Green Algal Models for Multicellularity.

The repeated evolution of multicellularity across the tree of life has profoundly affected the ecology and evolution of nearly all life on Earth. Many of these origins were in different groups of photosynthetic eukaryotes, or algae. Here, we review the evolution and genetics of multicellularity in several groups of green algae, which include the closest relatives of land plants. These include millimeter-scale, motile spheroids of up to 50,000 cells in the volvocine algae; decimeter-scale seaweeds in the genus Ulva (sea lettuce); and very plantlike, meter-scale freshwater algae in the genus Chara (stoneworts). We also describe algae in the genus Caulerpa, which are giant, multinucleate, morphologically complex single cells. In each case, we review the life cycle, phylogeny, and genetics of traits relevant to the evolution of multicellularity, and genetic and genomic resources available for the group in question. Finally, we suggest routes toward developing these groups as model organisms for the evolution of multicellularity.

Phototoxicities Caused by Continuous Light Exposure Were Not Induced in Retinal Ganglion Cells Transduced by an Optogenetic Gene.

The death of photoreceptor cells is induced by continuous light exposure. However, it is unclear whether light damage was induced in retinal ganglion cells with photosensitivity by transduction of optogenetic genes. In this study, we evaluated the phototoxicities of continuous light exposure on retinal ganglion cells after transduction of the optogenetic gene mVChR1 using an adeno-associated virus vector. Rats were exposed to continuous light for a week, and visually evoked potentials (VEPs) were recorded. The intensities of continuous light (500, 1000, 3000, and 5000 lx) increased substantially after VEP recordings. After the final recording of VEPs, retinal ganglion cells (RGCs) were retrogradely labeled with a fluorescein tracer, FluoroGold, and the number of retinal ganglion cells was counted under a fluorescent microscope. There was no significant reduction in the amplitudes of VEPs and the number of RGCs after exposure to any light intensity. These results indicated that RGCs were photosensitive after the transduction of optogenetic genes and did not induce any phototoxicity by continuous light exposure.

Three genomes in the algal genus Volvox reveal the fate of a haploid sex-determining region after a transition to homothallism.

Transitions between separate sexes (dioecy) and other mating systems are common across eukaryotes. Here, we study a change in a haploid dioecious green algal species with male- and female-determining chromosomes (U and V). The genus Volvox is an oogamous (with large, immotile female gametes and small, motile male gametes) and includes both heterothallic species (with distinct male and female genotypes, associated with a mating-type system that prevents fusion of gametes of the same sex) and homothallic species (bisexual, with the ability to self-fertilize). We date the origin of an expanded sex-determining region (SDR) in Volvox to at least 75 Mya, suggesting that homothallism represents a breakdown of dioecy (heterothallism). We investigated the involvement of the SDR of the U and V chromosomes in this transition. Using de novo whole-genome sequences, we identified a heteromorphic SDR of ca 1 Mbp in male and female genotypes of the heterothallic species Volvox reticuliferus and a homologous region (SDLR) in the closely related homothallic species Volvox africanus, which retained several different hallmark features of an SDR. The V. africanus SDLR includes a large region resembling the female SDR of the presumptive heterothallic ancestor, whereas most genes from the male SDR are absent. However, we found a multicopy array of the male-determining gene, MID, in a different genomic location from the SDLR. Thus, in V. africanus, an ancestrally female genotype may have acquired MID and thereby gained male traits.

Co-expression networks in Chlamydomonas reveal significant rhythmicity in batch cultures and empower gene function discovery.

The unicellular green alga Chlamydomonas reinhardtii is a choice reference system for the study of photosynthesis and chloroplast metabolism, cilium assembly and function, lipid and starch metabolism, and metal homeostasis. Despite decades of research, the functions of thousands of genes remain largely unknown, and new approaches are needed to categorically assign genes to cellular pathways. Growing collections of transcriptome and proteome data now allow a systematic approach based on integrative co-expression analysis. We used a dataset comprising 518 deep transcriptome samples derived from 58 independent experiments to identify potential co-expression relationships between genes. We visualized co-expression potential with the R package corrplot, to easily assess co-expression and anti-correlation between genes. We extracted several hundred high-confidence genes at the intersection of multiple curated lists involved in cilia, cell division, and photosynthesis, illustrating the power of our method. Surprisingly, Chlamydomonas experiments retained a significant rhythmic component across the transcriptome, suggesting an underappreciated variable during sample collection, even in samples collected in constant light. Our results therefore document substantial residual synchronization in batch cultures, contrary to assumptions of asynchrony. We provide step-by-step protocols for the analysis of co-expression across transcriptome data sets from Chlamydomonas and other species to help foster gene function discovery.

Volvox barberi (Chlorophyceae) actively forms two-dimensional flocks in culture.

Volvox barberi is a multicellular green alga forming spherical colonies of 10,000-50,000 differentiated somatic and germ cells. We observed that in culture, these colonies actively self-organized in just a few minutes into "flocks" that contained as many as 100 colonies moving and rotating collectively for hours. The colonies in flocks formed two-dimensional, irregular, active crystals, that is, geometric lattices within which individual colonies rotated separately. These groupings sometimes disassembled back into individual colonies just as quickly, but in some cases, flocks persisted over several hours. Close inspection of flock formation in the presence of a tracer dye suggested that colony and flock rotations were producing vortices in the fluid medium over a range spanning multiple flock diameters, perhaps providing a physical mechanism for aggregation.

Identification and characterisation of the Volvox carteri Moco carrier protein.

The molybdenum cofactor (Moco) is a redox active prosthetic group found in the active site of Moco-dependent enzymes (Mo-enzymes). As Moco and its intermediates are highly sensitive towards oxidative damage, these are believed to be permanently protein bound during synthesis and upon maturation. As a major component of the plant Moco transfer and storage system, proteins have been identified that are capable of Moco binding and release but do not possess Moco-dependent enzymatic activities. The first protein found to possess these properties was the Moco carrier protein (MCP) from the green alga Chlamydomonas reinhardtii. Here, we describe the identification and biochemical characterisation of the Volvox carteri (V. carteri) MCP and, for the first time, employ a comparative analysis to elucidate the principles behind MCP Moco binding. Doing so identified a sequence region of low homology amongst the existing MCPs, which we showed to be essential for Moco binding to V. carteri MCP.

Targeted migration of pherophorin-S indicates extensive extracellular matrix dynamics in Volvox carteri.

Hydroxyproline-rich glycoproteins (HRGPs) constitute a major group of proteins of the extracellular matrix (ECM). The multicellular green alga Volvox carteri is a suitable model organism in which to study the evolutionary transition to multicellularity, including the basic principles and characteristics of an ECM. In Volvox, the ECM is dominated by a single HRGP family: the pherophorins. Our inventory amounts to 117 pherophorin-related genes in V. carteri. We focused on a pherophorin with an unexpected characteristic: pherophorin-S is a soluble, non-cross-linked ECM protein. Using transformants expressing a YFP-tagged pherophorin-S we observed the synthesis and secretion of pherophorin-S by somatic cells in vivo, and we then traced the protein during its conspicuous migration to the ECM around prehatching juveniles and its localized concentration there. Our results provide insights into how an ECM zone surrounding the progeny is remotely affected by distantly located parental somatic cells. In view of the properties and migration of pherophorin-S, we conclude that pherophorin-S is likely to act as an ECM plasticizer to allow for dynamic ECM remodeling.